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Abstract 

This work is devoted to the formulation and derivation of the r]— /r/gamma and A—///gamma dis¬ 
tributions which correspond to physical fading models. These distributions are composite and are based 
on the 77 — p and A — // generalized multipath models, respectively, and the gamma shadowing model. 
Novel analytic expressions are derived for the corresponding envelope probability density functions. 
Importantly, the proposed models provide accurate characterisation of the simultaneous occurrence of 
multipath fading and shadowing effects which is achieved thanks to the remarkable flexibility offered by 
their parameters that render them capable of providing good fittings to experimental data associated with 
realistic communication scenarios. This is additionally justified by the fact that they include as special 
cases the widely known fading models such as Hoyt/gamma, Nakagami-m/gamma and Rayleigh/gamma. 
As a result, they can be meaningfully utilized in various analytical studies related to the performance 
evaluation of digital communications over composite multipath/shadowing fading channels. 

Index Terms 

p — p Distribution, A - /z distribution, gamma distribution. Multipath fading, shadowing, composite 
fading channels. 
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I. Introduction 

It is widely accepted that fading is a physical phenomenon that degrades communication 
signals during wireless propagation. A common approach for accounting for this effect is to 
exploit appropriate statistical distributions which are typically known as fading models. As 
a result, statistical models such as Rayleigh, Nakagami-m, Weibull and Nakagami-q (Hoyt) 
have been shown to be capable of modelling small-scale fading in Non-Line-of-Sight (NLOS) 
communication scenarios, whereas the Nakagami-n (Rice) distribution has been typically utilized 
in characterizing multipath fading in Line-of-Sight (LOS) communication scenarios, |fTfl-lf3l 
and the references therein. Capitalizing on the these models, M. D. Yacoub initially proposed 
three generalised fading distributions, namely, the a — /i, the n — //, the // — // models and 
subsequently the A — n and the k — // Extreme models, ffl- llTOll . These models are distinct thanks 
to the remarkable flexibility offered by their named parameters which render them capable of 
providing adequate fitting to experimental data which correspond to realistic communication 
scenarios. Their usefulness is also evident by the fact that they include as special cases all the 
aforementioned small-scale fading distributions. 

However, it is recalled that a crucial aspect of wireless radio propagation is that the occurrence 
of multipath and shadowing effects is typically simultaneous. As a consequence, in spite of the 
undoubted usefulness of the aforementioned fading models, they all ultimately fail to account 
concurrently for both shadowing and multipath fading. In other words, the utilization of the 
aforementioned fading models is limited to the effective characterisation of either the one or 
the other effect. Based on this principal limitation, the need for composite statistical models 
which are capable of providing efficient characterization of fading ”as a whole”, became evident 

m-m. 

To this effect, the authors in ifTTl proposed the Rayleigh/gamma fading model, which is also 
known as AT-distibution, (A'). Likewise, Shankar in lH2l exploited the flexibility of Nakagami-m 
distribution, which includes Rayleigh distribution as a special case, and introduced the Nakagami- 
m/gamma composite distribution - or generalised A'-distribution, (Kg). In the same context, 
the Weibull/gamma composite distribution was proposed in lfT3ll while an introduction to more 
generalized composite distributions was reported in fl4ll - llT9ll . 

Capitalizing on the above, the aim of this work is the formulation and derivation of the 
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77-77/gamma, and the A-77/gamma composite fading distributions. Specifically, after formulating 
these models, novel analytic expressions are derived for their corresponding envelope probability 
density function (pdf). The validity of the offered results is justified numerically and their 
behaviour is examined under different parametric scenarios. Thanks to their relatively convenient 
algebraic representation, the offered results are considered useful mathematical tools that can 
be efficiently utilized in analytical studies related to the performance of digital communications 
over 77 — ///gamma and A — ///gamma composite fading channels. To this effect, they can be 
meaningfully exploited in the derivation of analytic expressions for critical performance metrics 
such as error probability, channel capacity and higher order statistics, among others. 

The remainder of this paper is organised as follows: Section II revisits the basic principles 
of the 77 — 77, the A — 77 and the gamma distributions. Sections III and IV are devoted to 
the presentation, formulation and derivation of the 77 — 77/gamma and A — 7//gamma fading 
models, respectively, along with the necessary analysis on their behaviour. Finally, discussions 
on their potential applicability in wireless communications along with closing remarks are given 
in Section V. 


II. The 77 - 77, A - 77 and gamma Fading Distributions 


A. The 77 — 7/ Fading Model 

The 77 — 77 distribution is a general fading distribution that accounts for the small-scale variation 
of fading signals in NLOS communication scenarios. It is written in terms of two physical 
parameters, 77 and 77; the former corresponds to the ratio of the powers between the multipath 
waves in the in-phase and quadrature components, whereas the latter is related to the number 
of multipath clusters in the environment H, [HO- For a fading signal with envelope R and 
f = y/E(R 2 ), the 77 — 77 envelope probability density function (pdf) is mathematically expressed 
as 



( 1 ) 


where E{.) and r denote expectation and the root-mean-square ( rms ) value of R, respectively. 
Furthermore, the parameter 77 is positive and is given by 
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( 2 ) 


where var(.) denotes variance and the parameters h and H are defined as, 


2 + 7 / 1 + 7 / 


(3) 


4 


and 



(4) 


Also 0 < 7 / < oo is as already mentioned the scattered wave power ratio between the in-phase 
and quadrature components of each cluster of multipath. By taking the ratio of h and H, it is 
shown that [® 


Equation ([!]) can be also expressed in a normalized form by considering a normalized envelope 
P = R/f. Based on this, the corresponding normalized power pdf can be deduced by setting 

Pw(w) = p P (y/w)/ 2 y/w, namely, 



e- 2 '"”V i ( 2 nHw) 


( 6 ) 


It is recalled that the p — p model includes as special cases other widely known small- 
scale fading distributions. More specifically, the Hoyt distribution is obtained for // = 0.5 and 
b = (1 — 77)/(I + 77) (where b denotes the Nakagami-g or Hoyt parameter). Likewise, the one sided 
Gaussian distribution is attained for p —> 00 or 77 —>• 0 whereas the Nakagami-m distribution is 
attained for p = m and p —> 0 or p —> 00. Finally, the p — p distribution reduces to the Rayleigh 
distribution for the special case p = 0.5 and p = 1 , [ 0 . 

B. The A — // Fading Model 

The A —p distribution is a fading model which constitutes the Format -2 of the p—p distribution 
in HI. Recalling the same initial conditions as in the previous Subsection, its normalized envelope 
pdf is obtained by setting in (□}, 
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(7) 


and 



( 8 ) 


Therefore, for r] = A, it follows that 



(9) 


which is valid for — 1 < A < 1 . To this effect, the corresponding pdf is given by |j 9 ] eq. ( 1 )], 
namely, 



( 10 ) 


where 



(ID 


By utilizing once more the relationship pw(u>) = Pp(V™)/ 2 \/w, its normalized power pdf is 
expressed as 



( 12 ) 


Physically, likewise the rj—p model (Format 1 ), the parameter // is the inverse of the normalised 
variance and relates to the number of multipath clusters. On the contrary, the parameter A, or rj 
in Format 2 of the rj — // model, denotes the correlation coefficient between the scattered wave 
in-phase and quadrature components of each cluster of multipath (S). Importantly, the Nakagami- 
m fading distribution is obtained for /; = m and Ail, whereas for the special case // = 0.5 
and A = 1 , the A — n distribution reduces to the Rayleigh distribution 0 . 
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C. The Gamma Fading Model 

The log-normal distribution has been largely considered the optimum statistical model for 
characterising the shadowing effect, 111 - 0 . Nevertheless, in spite of its usefulness, it has been 
largely shown that when it becomes algebraically involved with other elementary and/or special 
function, its algebraic representation often renders it inconvenient to handle analytically. This is 
particularly the case in studies related to the analytical derivation of critical performance measures 
in digital communications over fading channels. Motivated by this, the authors in [fill proposed 
the gamma distribution as an accurate substitute to log-normal distribution. Mathematically, the 
envelope pdf of gamma distribution is given by [ 11 , eq. ( 4 )], namely, 


Py(v) = 


6-1 - M - 
y° L e n 


y> o 


( 13 ) 


r{b)n b ’ 

where the term b > 0 is its shaping parameter and Q = E{Y 2 ). This fading model has been shown 
to provide adequate fitting to experimental data that correspond to realistic fading conditions. 
In addition, its algebraic representation is particularly tractable and therefore, easy to handle 
both analytically and numerically. As a result, it has been widely considered a useful model for 
characterising shadowing effect and based on this, it has been exploited in the formulation of 
the K and K c composite multipath/shadowing models, [HTTl . Ifl 2 ll . 


III. The r] - ///gamma Fading Distribution 

A. Model Formulation 

According to the basic principles of statistics, the envelope pdf of a composite statistical 
distribution is constituted by superimposing two or more statistical distributions. In the present 
case, this is realized by superimposing one multipath and one shadowing distribution, namely, 

/*oo 

p R ( r )- p RlY ( r I y)p Y (y) d y ( 14 ) 

Jo 

where p RlY (x \ y) denotes the corresponding multipath distribution with mode y. Evidently, the 
rj — ///gamma composite fading distribution is formulated by firstly setting r — x and f — y in 
(|H) and then substituting in (fl4l) along with equation (fl3l) . To this end, it immediately follows 
that 
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= 4 f 00 e ( 2 ^ 2 ) 

x r(/x)^-lr(6)o ft Jo y 2 ^+i )-b e w£ 

Importantly, the term y 2 in (fl5l) has emerged from the term f 2 which denotes that the mean- 
squared value of the fading amplitude follows the gamma distribution. However, it is noted here 
that it can be also assumed that the root-mean-squared value of the fading amplitude is gamma 
distributed. In fact, this is exactly the difference between the Rayleigh/Lognormal and Suzuki 
fading models since in the former the runs value of the fading amplitude is modelled as by log¬ 
normal distribution, whereas in the latter it is the mean-squared value of this amplitude which 
is assumed to be log-normally distributed |J3). Therefore, by applying this principle in (fl5l) and 
letting u — 1/y, it follows that y = 1/u and du/dy = — 1/y 2 . To this effect and by taking the 
corresponding absolute value, the following expression is deduced, 


p x {x) 


4r°°e on / M _1 (2pHx 2 u) 

T(p)H^T(b)CL b Jo u b ~^ e 2fihx2u 


( 16 ) 


B. The Special Case p G N 

A closed-form expression for the envelope pdf of the // — ///gamma distribution for the case 
that /i 6N can be obtained with the aid of the closed-form series representation for the I n {x) 
function in [ED eq. (8.467)], namely, 


n 

f n+ i(a:) = ^ 


(n + k)\ [(-l) fc e* + (-l) n+1 e“*] 


k =0 


y/7rk\(n — k)\(2x 


\ k+ h 


n e N 


As a result, by performing the necessary change of variables, recalling that T(x) = (x 
substituting in (fT6l) yields, 


(17) 
1)! and 


T(p + k)p^ k h fJ 'x 2 (> 1 1 

Px( x ) = }L 


^ k\T(p - k)2 2k+1 r(p)r(b)H»+ k n b 


k =0 


f /'UU /'UU 

x \ (-l) k / u^- b - k - 1 e-^e~ 2flx ^ h - H)u du+ (-1)^ / w M-&-fc-i e -Tfe e -2/^ 2 (fc+ff)^ 


( 18 ) 

The above integrals can be solved in closed-form with the aid of fl2TL eq. (3.471.9)], namely, 








(19) 


J x v x e *e lx dx = 2 ^ K u (2a/ 07 

where K v (x) denotes the modified Bessel function of the second kind, l22l . Therefore, by 
performing the necessary variable transformation and substituting in (fl8l) . one obtains 


g+sk-fc B±%dL, a , n v-k+b-l 

Px (x) = S8(x) + y - g ^ )khX -— 

y o k\T(m - k)T(b)H»+ k n b+I ^- 

2 ,/spr) 


x 


I^/i—b—k I 2x 


2 fi(h-H) 


K„, 


fi—b—k 


(-1 )~ k (h-H)^ 


+ 


. . , , _ Li—b—k 

(—+ H) 2 


( 20 ) 


where 5(.) is the Dirac delta function and (x) n = Y(x+n)/Y(x) is the Pochhammer symbol f[2Tll . 
Furthermore, the parameter S denotes a normalisation scalar constant that must be determined 
so that equation (l20l) constitutes a true pdf. To this end, by recalling that / 0 °° p x (x)dx = 1 and 
/ 0 °° 8(x)dx = 1, one obtains, 


U—l t) ii+3k-b ft+b-k . . 

s = 1 _j2 2 2 v 2 


y k\T(m - k)Y{b)H»+ k Yl b+l ^ 


f-cxD 1 '-/i-b-k 


(h - H)^ Jo 


Kn-h-k I , { _ ir 

doc I 


foo Kfi—b—k 


2 l /sp>) 




1 +fc—/x —6 


(/l + if)*"^ J 0 


X 




( 21 ) 


Importantly, the above integrals can be solved with the aid of ll2Tl eq. (6.561.16)], namely, 


’ a ts n U 2 al fa + l+n\ fa + l-n 
x K n (bx)dx = —Y ( --- I r 


( 22 ) 


which is valid for Re b > 0 and Re a +1 ±n > 0. To this effect, by making the necessary change 
of variables and substituting in (1271) yields the following relationship, 


11-1 


s = i-E 


(/U)fc/l M 


^ k\{2HY+ k 


-If 


+ 


-1 y 


(23) 


( h-Hy- b ~ k (h + Hy- b ~ k _ 

Therefore, by substituting (l23l) into ([2Tb . a closed-form expression for the envelope pdf of the 
//-///gamma fading model for // e N is finally deduced, namely, 


dx 




































Px( x ) = 


M-1 1 | K^-b-k I 2x 


ftII ■ h—b I 2(T 


2 Kh-H) . 

- n - J l\^-b-k 

+ 


2fj,(h+H) 

n 


fc!T(m - k)T(b)H^+ k n b+!L ^ I (-l)-^-#)^ (-l)-/*(/i+ #)*=$"* 


k\{2HY+ k 


■ 1 )* 


m {h-HY 


—b—k 


+ 


-i y 


(.h + HY 


—b—k 


+ 1 


(24) 


C. 77ie general case fieK 

By recalling that the the root-mean-squared value - and not the mean-squared value- of the 
fading amplitude is gamma distributed, equation (fl5l) is alternatively expressed as, 


p x (x) = 


4^/n n^ + 2 h^x 2 ^ ft 

Y(fi)H^Y(b)n b Jo 


e~* Vj ( 2 ^y 

,i+u- bp ^ h f 


dy 


(25) 


Notably, the modified Bessel function can be also expressed in terms of the polynomial approx¬ 
imation in [23] eq. (19)] as follows: 


I ( X )~Y _ F(n + 0 _ — _f-i ^ (26) 

_ w r(f + l)r(n-(+l)r(i/ + f + l) '2/ ' ' ’ 

As n —* oo, the above expression reduces to the infinite series in [21, eq. (8.445)]. Therefore, 
by making the necessary variable transformation and substituting in (l26l) . one obtains 


Px( x ) = 


r (n + k)/u 2 ^ + ftH 2 k hfti 1 ~ 2 k x 4 ^ + ^- 1 n~ b 
—' r (k + l)T(n — k + l)r (/r + k + 2) T {fi)T(b) 


4y/n 


b— 0 ,, — 0 b — -\ — y. — 2 fthx 

y b ^ 1 e n e y 


k =0 


dy 

(27) 


By utilizing again (fl9l) in (1271) and after long algebraic manipulations, the following closed-form 
expression is deduced, 


Px( x ) 


k+3 y/nT(n + k)yX +tl+k H 2k 
* + £^ k\Y{n-k + l)Y{yb + k + \)Y(y) 


jjl —2/c^,2(/z-|-/c)-|-6— 1 

h k ~^Y(b)n^ +k Kh ~ 2{tx+k) 



(28) 


The scalar parameter S that constitutes p x (x) a true pdf can be determined by following the 
same methodology as in the previous Subsection. To this end, it is noticed that the corresponding 
resulted integral has the same algebraic form as the integral in ([22]) . Based on this, by making the 
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necessary change of variables, substituting in (l28l) and recalling that J 0 °° p{x)dx = J 0 °° 8(x)dx = 
1 yields the following expression 

_ v^r(n + k)H 2k n l ~ 2k T{2n + 2/c)2 1 - 2 ^ +fc ) 

fc!r(n-fc + i)r(/r + fc + |)r(/i)^+ 2fc ' 

As a result, by substituting the above expression into (1281) . one obtains the following analytic 
expression for /igM, 


Px( x ) 


n 

=i-E 


fc =0 


vAFr(n + k)H 2k n 1 ~ 2k Y{2fi + 2k)2~ 2 ^ +k 1 
2k\T(n -k + l)T(ii + k + l) r(/x)^+ 2fc 


22 ^ fc+3 y / 7rr(n + k)ii^ +fl+k H 2k j n 1 1 

+ fc!r(n-fc + i)r(/x + fc + i)r(/x) {/rMr(&)^+ fcA '" 2(At+fc) 



(30) 


To the best of the authors’ knowledge, expressions in (l24l) and (l30l) . have not been previously 
reported in the open literature. 


IV. The A - /i/gamma Fading Distribution 

A. Model Formulation 

As already mentioned in Section I, the A — n fading model is a different version of the rj — // 
fading model (format 2). As a result, the algebraic form of the corresponding envelope pdf of the 
two distributions is the same and therefore, analytic expressions for the A — /i/gamma composite 
model can be derived based in the same manner as the rj — /i/gamma model in Section III. To 
this end, by assuming that the rms of the fading amplitude follows the gamma distribution, 
setting r = x in ( ITOl) and substituting it in (fl4l) along with (fTTI) . yields an explicit expression 
for the envelope pdf of the A — /i/gamma distribution, namely 

= /■” e ~ Sj -4 G^ir) , (31) 

Px r(ii)VT=yr{b)n b Jo y i+„-b,,^T y 
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B. The special case / 16 N 

By representing the I n (x) function in (BTT) according to ( fT71) and setting u—l/y, one obtains 
the following expression, 


11-1 

p x ( x ) = 

k =0 


r(/i + k)^~ k (i - a 2 )** 2 ^-*)- 1 
k\T(p - k)2 2k + l T(p)T(b)\»+ k kl b 


x 


■ 1 )* 


u 


, 2ux J (l-A) 


du + (—1) A 


00 /II 1 2,j.x 2 (1+\) 

u n-b-k- lg-nsg T^x^ u du 


(32) 


l Jo Jo ) 

which has the same algebraic representation as (fl8l) . Therefore, by performing the necessary 

change of variables and making use of (fl9l) . the following expression for p x (x) is deduced, 


Px( x ) 


ii-i 

ss( X ) +yy 

k =0 


4^ - A2 )^ 

2 k\T{m — fc)r(6)A^ +fc f2 £ ^ 


x 


Kfi-b-k {^ X \J 

n b (-i)- fc (i- A)^ 


K,- b -k ( 2 xy^) | 

+ n b (-l)-^(l + A)^ J 


(33) 


Next, by following the same procedure as in Section III, the scalar parameter S is given by, 


M-l 


*=i-E 


) 2 - 


/^+3fc —6 /i,+6—fc , \ / 


u+fc — 6 

A 2 ) — 


/c =0 


fc!T(m - k)r(b)X^+ k n b+t 


X 


( —l) fc r°° K»-b-k {^ X \J 0(1 -A 2 )) ( — 1)^ f°° Kn-b-k {^ X \J fJ(i-A 2 )) 

( 1 _A)^ X ^ ^ + (l + A)^ Jo x i+k-»-b x 

(34 

Importantly, the above two integrals belong to the same class as the integral in (l22l) . Therefore, 
by performing the necessary variable transformation it follows straightforwardly that, 

^ (p) k (f - A 2 ) fc+1 

fy o k\(2\y+ k 

where A = /i — b. Evidently, by substituting ([331) into (l34l) yields a closed-form expression for 
the envelope pdf of the A — p fading model for // e N, namely, 


-i r 


(l-A) 


A—k 


+ 


- 1 )' 


(1 + A) 


A-k 


(35) 
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M- 1 

Px(^) = X] 


k =0 


4/i £± ^( / 2) fc (1 - A 2 )^ x^- fc+fe -i 

2 i£± ¥^fc!r(/2 - fc)r(6)A^+ fc n^ 


IS ( O r / 2/j(l-A) '\ r - /o T / 2/x(l+A) '\ 

j^n-b-k y^ x y n(i-A 2 )y ^n-b-k y ^(i-a 2 )y 

O b (—l)- fc (l - A)^ + (T'(-l)-^(l + A)^ 


^ Qjfc (1 - A 2 ) fc+1 
^ fc!(2A) M+fe 


(~l) fc (~1) M 

(1 - A)- 4 ^ (1 + A)^" fc 


(36) 



Fig. 1. Envelope pdf of the 77 -p/gamma distribution for 6 = 1.2, Q = 0.8, = 0.6 and different values of 77 . 


C. 77ie general case /i£l 

An explicit expression for the envelope pdf of the A — /i/gamma distribution when ft e E, can 
be derived with the aid of the polynomial approximation for the I n (x) in ( 1261) . To this effect, 
by making the necessary variable transformation and inserting in (fTil) yields 

T(n + \ 2 k n l- 2 k x ±(ii+k)-l 

Px{X - t^ 0 k\T(n - k + 1) (1 - A 2 Y +2k T(fj, + k + ±) 
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X 




b-2n-2k-l-% 


2 px* 


e n e y^-^dy 


(37) 


r(/i)r(&)n* J 0 

Evidently, the above relationship can be expressed in explicit form by evaluating the involved 
integral. Likewise in previous Sections, the algebraic representation of this integral is the same as 
the integral in (fl9l) . As a result, the following analytic expression is deduced in a straightforward 
manner, 


p x (x) = S5(x) + J2 


22 v k+3 y/nT(n + k)ii2 + ^ +k \ 2k 


n 


k =o 

l-2k x 2(n+k)+b-l 


k\T(n — k + l)r [p + k + |) T(p) 


x 


(i-A’r*r(&)nt+^“~^ V v 

The scalar normalization constant S needs to be also determined in order to constitute p x (x) a 
true pdf, i.e. / 0 °° p x (x)dx = 1. To this end, by integrating both parts of (l38l) from zero to infinity, 
recalling that J 0 °° 8{x)dx = 1 and utilizing (l22l) yields the following closed-form expression 


-K h _ 


b—2(fi-\-k) 



2p 


(38) 


s=i-E 


2 v/iT(n + k)n l - 2k \ 2k (1 - \ 2 y T( 2p + 2k) 


(39) 


t'o 22{fl+k)k '- T (n - k + 1)T (/i + k + i) T(p) 

Finally, by substituting ([39b in (l38l) yields an analytic relationship for the envelope pdf of the 
A — /i/gamma fading distribution for the case of p 6 R, namely, 


Px(x) = 1 + Y^ 


k=0 


1 — 2 fc^, 2 (/r+fc)+ 6—1 


2 |-M-fc+ 3 v /ET(n + k)p^ l+k X 2k 
—' k\T(b)T(n - k + 1)T [p + k + |) T(p) 


x 


-Kh- 


(i - a^)‘ + 4 VI «(i - v) 

V7rT(n + k)v}- 2k X 2k (1 - A 2 ) M T(2/j + 2k) 

~ ^ 2 2 (^+ fe )-ifc!r(n - k + i)r (p + k + ±) r(p) 



2 p 


To the best of the authors’ knowledge, the proposed expressions in equations 
novel. 


(40) 

and (l40l) are 














14 



Fig. 2. Envelope pdf of the 77 — /i/gamma distribution for b = 1.2, It = 0.8, 77 = 10 and different values of /r 

V. Numerical Results and Discussions 

In this Section, the general behaviour of the derived analytic expressions for the envelope pdf 
of the 77 — ///gamma and A — ///gamma fading distributions is demonstrated. To this end, Figure 1 
illustrates the pdf of 17 — /i/gamma with respect to x for b = 1.2, Q = 0.8, /./ = 0.6 and different 
values of 77. Likewise, Figure 2 considers b = 1.2, fl = 0.8, 77 = 0.6 and different values of 77 . 
In the same context, Figures 3 and 4 demonstrate the pdf of A — 77 /gamma for b — Q — 1.0, 
77 = 0.6 and different values of A and b = 1.25, = 1.5, A = 0.5 and different values of 77, 
respectively. One can observe the flexibility of the proposed models which render them capable 
of providing adequate fittings to experimental results. 

A. Usefulness and Applicability in Wireless Communications 

It is widely known that the algebraic representation of crucial performance measures is 
critical in studies related to analytical performance evaluation of digital communications. This is 
obvious by the fact that when the algebraic form of a corresponding measure is convenient, 
significantly increases the possibility that the derived relationships be expressed in closed- 
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Fig. 3. Envelope pdf of the 77 — /i/gamma distribution for b = ft = 1.0, fj, = 0.6 and different values of A 

form. Therefore, the fact that the form of the offered analytic result have a relatively simple 
representation, renders this model convenient to handle both analytically and numerically. To 
this effect, the derived expressions can be efficiently applied in various analytic studies relating 
to the performance evaluation of digital communications over composite multipath/shadowing 
fading channels. Indicatively, the offered expressions can be straightforwardly utilized in deriving 
explicit expressions for important performance measures such as, error probability, probability of 
outage, ergodic capacity, channel capacity under different constraints and higher order statistics. It 
is recalled here that expressions corresponding to the aforementioned measures can obviously be 
derived in both classical and emerging technologies such as SISO and MIMO systems, diversity 
systems, cognitive radio and cooperative systems and optical communications among others. 

VI. Closing Remarks 

This work was devoted in the introduction, formulation and derivation of the r) — ///gamma 
and A —/i/gamma fading distributions. These models are formulated from the superimpose of the 
77 — /x and A — /i generalised small-scale fading models, respectively, and the gamma shadowing 
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Fig. 4. Envelope pdf of the p — /i/gamma distribution for b = ft = 1.0, A = 0.5 and different values of p 


model. These distribution are particularly flexible and they include as special cases the widely 
known Hoyt, Nakagami-m and Rayleigh fading models. Novel analytic expressions were derived 
for the envelope probability density function which can be considered a useful mathematical 
tool in applications related analytical performance evaluation of digital communications over 
multipath/shadowing channels. 
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